Estimating second-order characteristics of 
inhomogeneous spatio-temporal point processes: 

influence of edge correction methods and intensity estimates 
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Abstract: We restrict our attention to space-time point pattern data for which we have 
a single realisation within a finite region. Second-order characteristics are used to analyse 
the spatio-temporal structure of the underlying point process. In particular, the space-time 
inhomogeneous pair correlation function and i\-function measure the spatio-temporal clus- 
tering/regularity and the spatio-temporal interaction, and thus help with model choice. Non- 
parametric estimators of the second-order characteristics require information from outside the 
study region, resulting to the so-called edge effects which have to be corrected, and depend on 
first-order characteristics, which have to be estimated in practice. Here, we extend classical 
edge correction factors to the spatio-temporal setting and compare the performance of the 
related estimators for stationary/non-stationary and/or isotropic/anisotropic point patterns. 
Then, we explore the influence of estimated intensity function on these estimators. 

Keywords: Edge correction; Intensity estimation; Point process; Reduced second mo- 
ment measure; Spatio-temporal data. 

1 Introduction 

Space-time point pattern data (also called spatio-t emporal po i nt pa ttern data) are increas- 
ingly available in a wide range of scientific settings ( Vere-Jonesl . I2009I ). Data-sets of this kind 
often consist of a single realisation of the underlying process. Usually, separate analyses the 
spatial and the temporal components are of limited value, because the scientific objectives 
of the analysis are to understand and to model the underlying spatio-temporally interacting 
stochasti c mechanisms . Generic methods f o r ana lysing such processes are grow ing; see for 
example iDiggld (J20061 ) , iDiggle and Gabriell (|20ld ) and ICressie and Wikld (J201ll ) . There is 
already an extensive literat ure on the us e of p oint process models in the specific field of 
seismology; see, for example, IZhuang et all (|2002l ) and reference s therein. There a re basically 
two ways for modeling space-time point patterns (JDiggld . 120061 : IVere-Joned . l2009i ) . The first 
is descriptive and aims at providing an empirical description of the data, especially from 
second-order characteristics. The second is mechanistic and aims at constructing parametric 
point process models by specifying parametric models for the conditional intensity function. 
Here, we will consider the former and analys es will be base d on extensions of the pair corre- 
lation function and the Ripley's if-function (JRipleyJ, 119771 ) to summarize a spatio-temporal 
point pattern and test hypotheses about it. 



E. Gabriel - Estimating second- order characteristics of spatio-temporal point processes 



In many applications the point patterns cannot be considered homogeneous. This is par- 
ticularly so in epidemiological studies where the observed point pattern is spatially and tempo- 
rally inhomogeneous, as the pattern of incidence of the disease reflects bo th the spatial distri- 
butio n of the population at risk and systematic temporal variation in risk (JGabriel and Diggld . 
20091). In the sp atial case, there are different approaches to deal with the inhomogeneity 



(jUlian et all . 120081 ) : some consist in considering subplot s of the study region and employin g 
means of the methods for homogeneous point processes ( Allard et all . J2Q0ll : iBrix et all . l200ll ). 
others have their own unde r lying theory, mainly based on the inhomogeneous iT-function 
proposed bv iBaddelev et all ( 2000 ) . The latter have been f or ins tance used in environmen- 
tal e pidemiology (JDiggle et all . 120071 ) . ecology ( Law et all . |2009j) or economy ( Arbia et all . 
20121). The inhomogeneous if -function has been extended to the spatio-temporal setting 
bv lGabriel and Diggld (2009j). Second-order characteristics are thus analysed from the spatio- 
temporal inhomogeneous if -function (STIK-function) or equivalently from the spatio-temporal 
pair c orrelation function under the assumption of second-o rder intensity re- weighted station- 
arity (JDiggle and Gabriell . 120101 ; iGabriel et al . 20101. 2013). Spatio -temporal separability of 
the STIK-function has been studied in lMoller and Ghorbanil (120121 ) . 

Here, we restrict our attention to space-time point processes for which we have a single 
realisation within a finite (spatio-temporal) region. Estimating the STIK-function or the 
pair correlation function thus faces two issues. On the one hand, they depend on first-order 
characteristics (see Section |2|) which are unknown in pr actice, but replacing the intensity by 
an estimate must be made carefully as it may imply bias ( Baddelev et all . l2000l ). The problems 
encountered when using the same point pattern to estimate both a spatially varying intensity 
and second-order characteristics can be overcome by adjus ting the estimate o f intensity so as 
to take account of exp lanatory variab l es; see for instance iDiggle et all ( 20071 ) in the context 
of case-control data or lMrkvicka et all ( 20121 ) for the est imation of c l uster ed poi nt processes 
with inhomogeneous cluster centres. As suggested in lArbia et all ( 20121 ) and IDiggle et al 
(J2007I ). o ther ways to get around this difficulty is to as sume a parametric model for the 
intensity (JMoller and Waagepetersenl . 120031 ; iLiu et all . 120071 ) or to allow separation of first- and 
second-order effects. The latter is particularly difficult to check in the absence of independent 
replications, as we cannot in general make an empirical distinction between first-order and 
second-order effects without additional assumptions. On the other hand, estimating the 
STIK-function or the pair correlation function requires information from outside the study 
region, resulting to the so-called edge-effects which have to be corrected. In the spatial case, 
the four main edge correction methods are: the border method (all sample points that are 
closer to the border than to their nearest neighbours are eliminated), the guard method (the 
sample points within the guard zone are not used for the analysis, but used as some nearest 
neighbours), the toroidal method (it makes eight identical copies of the rectangular study 
region and places them around the original one), the weighted method (or the isotropic 
method), also called Ripley's method (each weight is equal to the probability that points 
around the location of sample point i will be in the study region). There is an extensive 
literatur e on edge correc t ion for i nter-pqints di s tance methods in the spatial case; see fo r 



instance 



Baddelevl dl99flh 



Li and Zhangj tooit) a nd 




Diggld toah . llllian et all ifcood j. lLaw et all (J2009I ) 



for the definitions of usual edge correctio n factors 



Goreaud and Pelissierl (1199911. lHaasd (I1995T) and 
other specific ones, and lYamada and Rogersonl ( 



some 



Pommerening and Stoyanl (200a) for 
2003J) for an empirical comparat ive study 



in the homogeneous case. In the spatio-temporal case, ICronie and Sarkkal ( 201ll ) propose 
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some edge corre ction met h ods w hen considering a given parametric model for the point 
pattern. Finally. iBaddelevi ( 19991 ) points out that edge effects are more difficult to overcome 



in dimension greater than two. 

As above mentioned, estimating first- and second-order characteristics from a single re- 
alisation of point process and correcting edge effects have been widely studied in the spatial 
case. The aim of this paper is to discuss and explore the influence of these two issues in 
the spatio-temporal setting, particularly on non-parametric estimates of the STIK-function 
and the pair correlation function. The structure of the remainder of the article is as follows. 
Section [2] reviews some existing results concerning the second-order characteristics of inho- 
mogeneous spatio-temporal point processes and their estimation. In Section [3j we extend 
classical edge correction factors to the spatio-temporal setting and provide an empirical com- 
parative study of the related estimators. Section [4] discusses the influence of using an estimate 
of the intensity on the performance of the estimators from a simulation study. Section [5] is a 
short discussion. 

2 Description of second-order characteristics 

2.1 Spatio-temporal point processes 

Following iDiggle and Gabriell (j2010l ) and iGabriel and Diggld (J2009i ) , let us consider spatio- 



temporal point processes whose events are defined as countable sets of points x = (s, t) where 
s denotes location and t denotes time. In practice we observe n events Xi = (si,ti) within 
a bounded spatio-temporal region S x T C R 2 x R. Here, we focus on the case where 
additional variables describing characteristics of the events (marks) are not available. In 
the following, N(A) denotes the number of events in an arbitrary region A. For for mal 
definitions of point process chara cteristics, see for example iDalev and Vere-Jonesl (J2003I ) or 
Moller and Waagepetersenl (120031 ). 



2.2 First-order and second-order characteristics 

First-order characteristics are described by the intensity of the process, 

w , w , ,. E[N(dsx dt)] 

X(x) = X(s,t) = Inn ,,,,,, n , 

|ds|->0,|dt|-M) |os||ot| 

where ds defines a small spatial region around the location s, \ds\ is its area, dt is a small 
interval containing the time t, \dt\ is its length. If the intensity is constant, X(s,t) = X for all 
(s,t), then the process is called homogeneous or first-order stationary. 

The relationship between numbers of events in pairs of subregions within S x T is described 
by the second-order characteristics. The second-order intensity is defined as 

, 2(1>l0= lim mmm, 

\A\,\A'\^0 \A\\A'\ 

where A = ds x dt and A' = ds' x dt' are small cylinders containing the points x and x' 

respectively, x ^ x'. 

Equivalent descriptors of second-order characteristics include the covariance density, 

7(x,x') = \<2(x,x') — X(x)X(x') 
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and the rad ial distribution function or the point-pair correlation function ( Cressid . Il993l ; 
DiggleLl2003 ): 

9[,) X(x)X(x')' [) 

For a spatio-temporal Poisson process, the covariance density is identically zero and the pair 
correlation function is identically 1 . Larger or smaller values than these benchmarks therefore 
indicate informally how much more or less likely it is that a pair of events will occur at the 
specified locations than in a Poisson process with the same intensity. 

In practice, it is often difficult to estimate these moments and one may need some relaxing 
assumptions. 

2.3 Spatio-temporal stationarity 

A point process is first-order and second-order stationary in both space and time, if: 

X(s,t) = X and A 2 ((s, t), (s' , t')) = X 2 (s - s', t - t 1 ). 

A stationary spatio-temporal point process is isotropic if X 2 ((s,t), (V,t')) = X 2 (u,v), where 
u = \\s — s'\\ and v = \t — t'\ denote spatial and temporal distances respectively. 
Baddelev et all ( 200CJ ) introduce a weaker assumption in the purely s patial case, called second- 
order intensity reweighted stationarity. iGabriel and Diggld ( 20091 ) extend it to the spatio- 



temporal setting as follows: an isotropic spatio-temporal point process is second-order in- 
tensity reweighted stationary if its its pair correlation function depends only on the spatio- 
temporal difference vector (u,v). 

2.4 Separability 

A spatio-temporal point process is first-order separable if its intensity X(s, t) can be factorised 
as 

X(s,t) = X 1 {s)X 2 (t), for all (s,t) G S x T, 

where Ai and X 2 are non-negative functions. 

The pair correlation function is separable in space and time if 

g{{s,t),{ S ',t'))= gi { S , S ')g 2 (t,t'), 

where gi and g 2 are non-negative functions. Then for an isotropic second-order intensity 
reweighted stationary point process, we have 

g{{s,t),{s',t'))= gi (u)g 2 (v). 

Some diagnostic proc edures for checking hypothes es of second-order spatio-temporal separa- 
bility are proposed in lMeSller and Ghorbanil ( 20121 ). 



2.5 Spatio-temporal inhomogeneous K- function 

Second-order characteristics can also be described through an extension of Ripley's K- 
function . For a second-order inten sity reweighted stationary, isotropic, spatio-temporal point 
process. IGabriel and Diggld ( 20091 ) define the space-time inhomogeneous K- function by 



K(u,v) = 2ttI / g(u',v') u'du'dv'. (2) 

-v JO 
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For an inhomogeneous spatio-temporal Poisson process with intensity X(x), the second- 



order intensity is A2 



x,x 



X(x)X(x'), hence the pair correlation function g(u,v) is identi- 



cally 1 and the STIK-function is K(u, v) = 2iru 2 v. Thus, K(u, v) — 2nu 2 v can be used as 
a measure of the spatio-temporal aggregation or regularity, using an inhomogeneous Poisson 
process as a benchmark. Positive values of K(u, v) — 2iru 2 v indicate clustering, or aggre- 
gation, at spatial and temporal separations less than u and v, respectively, while negative 
values indicates regularity. 

The STIK-function can also be used as a measure of spatio-temporal interaction. Its 
separability into purely spatial and temporal component s indicates absence of intera ction. 
In this case the ratio K (u, v)/{Ks(u)Kt(v)} is constant (jMoller and Ghorbanil . l2012l ). It is 
identically 1 for a Poisson process. 



2.6 Estimation 

An unbiased estimator of the STIK-function, proposed in iGabriel and Diggld ( 2009 ) under 
the assumption of isotropy, is 



n 

^— ' ^— ' 111: 



1=1 jy^i 



Wij X(xi)X(xj 



-{||s»— «j||<u ; \ti-tj\<v}> 



(3) 



where Wij is an edge correction factor to deal with spatial-temporal edge effects. The proof 

of unbiasedness is given in Appendix 1 for the edge correction factors defined in Section 13.11 

The space-time pair correlation function defined in Equation (pQ) can be estimated by 



i=i j^i j 



t (v - \ti -tj\) 



X(xi)X(xj) 



(4) 



where Wij is defined in Equation ([3]) and k s (-), kt(-) are kernel functions with bandwidths h s 
and hj. Experien ce with pair correlation function estimation recommends box kernels, see 
Illian et all ifeood ). 



These estimators depend on an edge correction factor and assume that the intensity is 
known. In practice the intensity function is unknown and have to be estimated. So, the 
questions are: which method can be used to correct edge effects and can be used to estimate 
the intensity function? what are their influence on the performance of the STIK-function 
and the pair correlation function? This is explored in the following sections. 



3 Influence of edge correction methods 

The literature o n edge correction methods is mo re ex tensive in the spati a l case , see e.g. 
Baddelevl (| 19981 ). Illlian et all (I2008T). iRiplevl (1198811 and IStovan and Stovanl (Il994l\ than in 



highe r dimensions (jBaddelev et all . 
20ld ). 



1993 



Cronie and Sarkkal . l201ll ; lDiggle et al 



199ll ; |Jafari-Mamaghani et all . 



In this section we extend three classical spatial edge correction factors to the spatio- 
temporal setting and compare the performance of the related estimators of the second-order 
characteristics for stationary/non-stationary and/or isotropic/anisotropic point patterns. 
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3.1 Edge correction methods 

The simple st approach t o dea l with spatio-temporal edge effects consists in correcting them 
separately (JDiggle et all . Il995l ) . Thus, the edge correction factor in Equations ([3]) and (jl]) is 
the product of a spatial edge correction factor and a tem poral edge correc tion f actor. Here, w e 
consider spatial edge correction factors (as described in iBaddeleyl (J1999J ) and iRipleyl (|1988l )) 
running for any shape of S. This eliminates in particular the toroidal method. 

In the study of Section [32]we also consider the case where no edge correction is performed: 



?c 



u 



\SxT\ 



3.1.1 Isotropic edge correction method 



The weight is proportional to the product between the Ripley edge correction factor ([Riplevl . 
19771 ) and its analog in one-dimension: 



p„ Tl (*) ( S ) 

Wij = ISxTlwl/wlj', 



where the temporal edge correction factor w, 



(t) 



2\U 



(t) 



u 



1 if both ends of the interval of length 



tj\ centred at L lie within T and w,„ 



1/2 otherwise ( Diggle et all . Il995l ) and w^/ 



the proportion of the circumference of a circle centred at the location S{ with radius ||sj 
lying in S. 



i.i 



3.1.2 Border and modified border methods 

These methods restricts att ention to those events lying more than u units away from the 
boundary of S ( Diggld . Il979l ) and more than v units away from the boundary of T. Thus, for 
d(si,S) denoting the distance between Si and the boundary of S and d(ti,T) denoting the 
distance between ti and the boundary of T, we have for the border method 



W; 



Z)j=l 1 {d(s j ,S)>« ; d(tj,T)>v}/H x j) 

1 {d(s i ,S)>u ; d(U,T)>v} 



and 



w 



Pen I x \m 



ev\ 



i.i 



l{d(si,S)>u ; d(U,T)>v} 

for the modified border method, where Sq u and T e „ are the eroded spatial and temporal 
region respectively, obtained by trimming off a margin of width u and v from the border of 
the original region (jBaddelev and Turner! . |2000| ) . 



3.1.3 Translation method 

The weight i s the proportion of translations of (xi,Xj) which have both xi and Xj inside S xT 



: sec 



Qhserl ( 19831 ) for the spatial 



case: 



Win 



\sns s 



x iTnr*.. 



where S q 



Jl] |(_> I I t->Si— Sj I • I -"• 1 I J -ti—tj\l 

and T^-t- are the translated spatial and temporal regions. 
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3.2 Simulations 

We examine the perform ance of the STIK-f unction and the pair correlation function via 
Monte Carlo simulations. iGabriel et all ( 20131 ) covers many of the models encountered in ap- 



plications of point process methods to the study of spatio-temporal phenomena. In particular 
an inhomogeneous Poisson process is defined by: 

1. The number N(S x T) of events within the region S xT follows a Poisson distribution 
with mean f s J T X(s, t) dt ds. 

2. Given N(SxT) = n, the n events in S x T form an independent random sample from the 
distribution on S x T with probability density function f(s,t) = X(s,t)/J s J T X(s',t') dt'ds'. 

and a Poisson cluster process is defined by: 

1. Parents form a Poisson process with intensity v(x). 

2. The number of offspring per parent is a random variable N c with mean m c , realised 
independently for each parent. 

3. The positions and times of the offspring relative to their parents are independently 
and identically distributed according to a trivariate probability density function </?(•) 
on R 2 x R. 

4. The final process is composed of the superposition of the offspring only. 



Reali sations of these point processes can b e generated using the R (jR Development Core Teaml . 
2012T ) package stpp (JGabriel et all2013h . 



For the simulation study, we set S x T = [0, l] 3 and simulate processes with expected 
number of points E[iV(S x T)] = n = 375. We generate N s i m = 1000 realisations of 

(i) HPP(X): an homogeneous Poisson process with intensity A = 375, 

(ii) IPP(X(x);f3): inhomogeneous Poisson processes with intensity function 

n/3 3 
Xi(x) = Ai(s,t) = _ jWfZ e -Pi) ex P(Pi(sx + s y - t)), (5) 

with s = (s x , s y ) and j3\ G {1, 2}, or 

. , , n (1.25 + cos(/3 2S:r + 0.25)) (1.25 + cos(/3 2 t + 0.25)) 

A2(x) = 71— (6) 

(1.25 + ^ (sin(/3 2 + 0.25) - sin(/3 2 )) 

with/3 2 e {3,5}, 

(iii) PCP\{y, a, a, m c ): stationary Poisson cluster processes with intensity of parents v = 25. 
The spatial distribution of the offspring is a zero-mean bivariate isotropic normal dis- 

(2) 

tribution <f> 2 with standard deviation a S {0.025,0.05,0.1,0.15,0.2} and the temporal 
distribution is exponential with rate a = 0.2 and denoted £ a (t). The expected number 
of offspring per parent follows a Poisson distribution with mean m c = 15. 



E. Gabriel - Estimating second- order characteristics of spatio-temporal point processes 



This Poi sson cluster process is an interpretation of a spatio-temporal shot-noise Cox 
process (JMgiller and Diaz-Avalosl . bold ; llVfollerl . 12003 ) with residual process R(x) - 



- Xlrie* ( P( X ~ v)i wnere ^ is a stationary Poisson process in M 2 x K with intensity v 

(2) 
and (p is the density function ip(u,v) = 4> 2 (\\ u \\)£a(v). For such a process, we have 

g(u,v) = 1 + -ip*(p{u, v), where * denotes convolution and ip(u,v) = ip(—u,—v). Thus, 



a l llull 2 



and 



5(M) = 1 + 8^ eXP "V " OT 



K(u,v) = 2ttu 2 v H (1 — exp(— av)) I 1 — exp I ^ I 

2v V V 4tT z / 



(iv) PCP2{v, <r,a, m c , (, 2 ,9, lu): geometric anisotropic Poisson cluster processes, as defined 
in (hi) but with g(u, v) = go ( VuS _1 n f , v ) , where u G M. 2 is a row vector with transpose 
the function go is such that g is locally integrable (see lM0ller and Toftagerl . |2012j), 



0.' . 



S is a 2 x 2 symmetric positive definite matrix of the form S = lj Ugdiag(l, £ )Ug 

with Ue = ( . ,„/ A, ). Here 6 = 7r/4, w = 4, and the anisotropy factor is 

\ sm(6') cos(6») J 

C G {0.25, 0.5, v 7 ^, v / 0?75, l}. 

(v) PCPs(X(x), (3, a, a, m c ): non-stationary Poisson cluster processes, with intensity of par- 
ents defined by Equation ([5]) or ((6|) and clusters as defined in (iii). 

FigureQ]shows static display of a realisation of the inhomogeneous Poisson process IPP(Xi(x), 
j3\ = 2) (top left) and IPP{\2{x), $2 = 5) (middle top), the stationary Poisson cluster pro- 
cess with a = 0.025 (top right), the non-stationary Poisson cluster process PCP^,{\2{x),j32 = 
5,<T = 0.1) (bottom left) and the anisotropic Poisson cluster process with a = 0.1 and 
£ 2 = 0.0625 (middle bottom) and £ 2 = 0.5 (bottom right). The time is treated as a quantita- 
tive mark attached to each location, and the locations are plotted with the size and colour of 
the plotting symbol determined by the value of the mark. Dark and large dots correspond to 
recent events in time. F urther forms o f stati c and dynamic display of spatio-temporal point 
patterns can be found in lGabriel et all ( 20131 ). 



3.3 Results 

The pair correlation function estimation is performed on the N s i m = 1000 realised data 
sets of each point process. The STIK- function is only estimated for the processes (i)-(iii) 
and (v). In both cases we use the different edge correction methods; we denote by Kj and 
57, Kb and 'gs, Kmb and guB, Kt and gr, the estimators of the STIK-function and of 
the pair correlation function when using the isotropic, the border, the modified border, the 
translation edge correction factor respectively and K^ and g^ the estimators in the case of 
no edge correction. In this section, we assume that the intensity is known to put ahead the 
influence of edge correction methods. 
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Figure 1: Realisations of the inhomogeneous Poisson ct Poisson cluster processes. Dark and 
large dots correspond to recent events in time. 



3.4 Performance of the estimators relative to edge correction methods 

Performance is first measured through the empirical bias and mean squared error (MSE) 
which are computed for each spatial and temporal distances u and v in the sequence starting 
from 0.01 to 0.25 by increment of 0.01. 

All estimators of the STIK-function have very small empirical bias and MSE at small 
spatial and temporal distances (less than 0.1) and their differences are observed for larger 
distances. Figure [2] shows the MSE of the STIK-function estimated when using the isotropic 
(I), border (B), modified-border (MB), translation (T), none (N) edge correction methods. 
Spatial distances u are in abscissa, while temporal distances v are in grey: the darker, the 
greater, with the same range of values than u. For the Poisson processes (both homogeneous 
and inhomogeneous), Kt, Kmb and Kjy have negative bias, while Kj has a positive one. 
The bias of Kb is almost constant. The MSE of Kt has larger values than the others for 
large values of u and v. This is illustrated for IPP{\\ (x), fi\ = 1) on the first row of Figure[2j 
Nevertheless, in all cases bias and MSE have negligible values. Values of the bias and MSE 
slightly increase for the stationary Poisson cluster processes PCP\{<j) and particularly when 
the range of spatial clustering, a, is small. For such processes, the bias of Kb increases with u 
and v, Kj and Kt perform less well and better, respectively, than the others. The second row 
of Figure [2] shows the MSE of all estimators for the process PCP\{a = 0.05). For the non- 
stationary Poisson cluster processes PCPs(Xi(x), /3i,a), bias and MSE of all estimators are 
similar to the stationary case. For PCPs(X2(x),/32,cr), all estimators have positive bias. Bias 
and MSE increase with a and particularly for Kj and Kmb ■ The last two rows of Figure [2] 
illustrate the MSE of the estimators for the processes PCP3 (Ai (x) , (3\ = l,a = 0.05) (third 
row) and PCP 3 (\ 2 (x), (3 2 = 3, a = 0.05) (last row). 

All estimators of the pair correlation function have small empirical bias and MSE with the 
exception of the case of very small values of u as the estimators contain the term "l/-u". Bias 
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Figure 2: Mean squared error of the STIK-function from IPP(Xi(x), f3\ = 1) (first row), 
PCPi(a = 0.05) (second row), PCP 3 (Ai (x) , fit = l,a = 0.05) (third row) and 
PCP3(\2(x),f32 = 3, er = 0.05) (fourth row). Spatial distances u are in abscissa; 
temporal distances are in grey: the darker, the greater, with the same range of 
values than u. 



at very small spatial distances could be reduced by the reflection technique ( Doguwal . ll990l ). 
For the Poisson processes HPP and PCPi(X(x),/3), we get similar comments than from the 
STIK-function, except that bias and MSE of gMB are almost constant with u and that gj has 
a large variability at small distances, which increase with the inhomogeneity of X(x). For the 
Poisson cluster processes PC Pi (a), PCP2(cr,( 2 ) an d PCPs(X(x), f3, a), smaller the spatial 
dispersion of offspring a and stronger the anisotropy (i.e. greater £ 2 ), greater the bias and 
the MSE at small u. This is particularly so for gi, gs and gMB- The MSE of all estimators of 
the pair correlation function are illustrated in Figure [6] (Appendix 2) for the same processes 
as in Figure E] and for PCP 2 {a = 0.2, C 2 = 0.0625) and PCP 2 (a = 0.2, C 2 = 0.5). 

The relative performance between the STIK-function and the pair correlation function 
does not only depend on the point process; it also depends on the edge correction method used 
to estimated the second-order characteristics. We compare the performance of K and g from 
the absolute relative error (ARE), \T(u,v) —T(u,v)\/T(u,v), evaluated for all (u, v) and for 
each edge correction factor, with T = K or g. For the Poisson processes, the STIK-function 
perform better than the pair correlation function whatever the edge correction factor, except 
for the border and modified-border methods in the case of IPP(\2(x), f3 2 = 5) for which 
we get comparable performances of the two estimators. For the Poisson cluster processes 
PCP\{a) with a < 0.15, gi and gs have smaller ARE than Kj and Kb for all distances, 
but not for the others at large distances. When a > 0.15, this phenomenon is stronger 
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and Kt can globally be better than gx- We get similar results for PCP3(\i(x),(3i,a). For 
PCP3(\2(x), 02, cr), g have to be preferred whatever the edge correction method used in the 
estimation. 

To avoid problems associated with multiple testing at different spatial and temporal scales, 
we also compute integral deviation measures 



D 



(T(u, v) — T(u, v)) 2 du dv, 



(7) 



where the summary statistic T is either K or g. We then compute the relative efficiency 
of the overall D defined by 100 x min^ 14/14, where 14 is the empirical variance 14 = 

Djfc — Dk) , k £ {I, B, MB,T,N}, with Di\, the deviation measure evaluated 



1 



N„ 



E-'* sim 
i=l 



from the estimation of T for the ith simulation and the /cth edge correction method and 



Di 



*JV 



A', 



^-\jim D ik ^ Tables H] and [5] (see Appendix 2) give the relative efficiency of the 
overall D defined from the STIK-function and the pair correlation function respectively. They 
are summarized in Table Q] which provides, for each point process, the best edge correction 
method, i.e. the one for which the relative efficiency is 100. It shows that the border method 
is the most efficient for Poisson processes (all cases for K and homogeneous or with weak 
inhomogeneity for g) and for the non-stationary Poisson cluster processes PCP^^ix),^)- 
Then for inhomogeneous and/or clustered isotropic and anisotropic processes, the translation 
method is the most efficient. 





K 


9 


HPP{\) 


Border 


Border 


IPP{X{x),p) 


Border 


Translation 


PCPtia) 


Translation 


Translation 


pcp 2 {cjX 2 ) 


- 


Translation 


PCP 3 {\i(x),/3i=l,a) 


Translation 


Translation 


PCP 3 (X 1 (x),/3i=2,cr) 


Modified Border 


Translation 


PCP 3 (A 2 (x),/3 2 ,a) 


Border 


Border 



Table 1: Most efficient edge correction method according to the point process. 



3.5 Clustering detection relative to edge correction methods 

Let us now look at the ability of detecting clustering relative to the edge correction method. 
The pair correlation function is evaluated from the realisations of the Poisson cluster pro- 
cesses PCP\{a) and PCP3(\(x), f3,a). We say that there is a clustering tendency when g 
is greater than the upper envelope of the pair correlation function evaluated under the null 
hypothesis (absence of clustering), g up (u,v) = maxj^fc(n, v), i.e. evaluated for the Pois- 
son processes HPP(X) and IPP(\(x), (3). Thus, for all distances u and v and for each 
edge correction method (subscript k), we compute the probability of detecting clustering: 

The envelopes are built fronr 1000 sinrulations, thus 



Pk(u,v) 



A s 



E 



N. 



im -1 ^_^ 

1 - L {9ifc(«,^)>9up(«:'")} 



leading to reasonable values for the type I error probabilities (between 0.0 2 and 0.04 for K and 
betwe en 0. 04 and 0.1 for q accordi ng to the edge correction method); see lLoosmore and Ford 
(|2006l ) and iGrabarnik et all (|201ll ) for a guidance on the use of envelop tests and deviation 
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tests. The probability of detecting clustering is expected to be maximum for spatial dis- 
tances u less than ~ 2a corresponding to the size of the clusters defined by our Poisson 
cluster processes. Because of the range of spatial distances used to estimate the pair correla- 
tion (between 0.01 and 0.25), we restrict our analysis to a G {0.025,0.05,0.1}. Results from 
the STIK- function are omitted for conciseness. 

For the stationary Poisson cluster processes PCP\(a), the probability is maximum for 
spatial distances up to 0.07, 0.12 and 0.18 when a = 0.025, 0.05 and 0.1 respectively. This is 
observed for the border, modified-border and translation methods and when no edge correc- 
tion is done. The isotropic edge correction factor leads to less powerful results as a increases. 
We get similar results for PCP^{\\{x),j3i = 1, er), than for PCP\{a = 0.025). However, for 
9MB, 9t an d <?jv the probability of detecting clustering decreases both when a and v increase. 
This is also the case for PCP%(\\{x),f5\ = 2, a), but the probabilities are weaker for all u 
and v because, whatever the edge correction method used, g is under-smoothed for this pro- 
cess. This can be explained by some overlapped clusters. For the processes PCP3(\2,/32,cr), 
5/ is only able to detect clustering when a = 0.025. The estimators gs, gMB-, 9t and g^ 
behave similarly than for PCP^{X\{x),j3i = 1, a), except that "false detections" appear for 
large values of u. This is straightly related to over-smoothed second-order characteristics 
for such processes. To illustrate all these comments, the probability of detecting clustering 
of gs and g~T are plotted in Figure [3] for the Poisson cluster processes PCP\{a = 0.025), 
PCP\{(J = 0.05) and PCPs(\(x), /3, a) with a spatial dispersion a = 0.025. Grey shading 
indicates values of the probability of detecting clustering. Black corresponds to one and white 
corresponds to 0. 
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Figure 3: Probabilities of detecting clustering using the border method (first and third 
columns) and translation method (second and fourth columns) for the Poisson clus- 
ter processes PCPi(a) and PC P^(X(x) , (3 7 a) . 
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Similarly, we obtain probabilities of detecting clustering for the overall statistic D. How- 
ever, to avoid misleading results due to bias in the estimations, we replace D defined in 
Equation ([7]) by 



Afc 



o Jo 



(dik(u,v) -g k (u,v)) dudv, 



(8) 



for k £ {I,B,MB,T,N}, i = 1, . . . ,N sim and g k (u,v) 



^n. 



^ sir, 



Ei=1 m 9ik(u,v). Table [2] shows 



these probabilities for various scales of spatial clustering and inhomogeneity. The performance 
of the estimators are good for the stationary Poisson cluster processes PCP\{a). As soon 
as we introduce spatial inhomogeneity, gj is no longer able to detect clustering. The border 
method appear to provide the best results. For PCPs(X(x),/3, a), the power of each estimator 
is related to the range of spatial clustering a, it decreases as a increases, and to the model 
of inhomogeneity X(x). 







Edge correction factor 




Point process 




I B 


MB 


T 


N 


PCP 1 {a = 0.025) 




1.000 1.000 


1.000 


1.000 


1.000 


PC7Pi(<t = 0.05) 




0.856 1.000 


1.000 


1.000 


1.000 


PCPi(cr = 0.1) 




0.007 1.000 


0.999 


1.000 


1.000 


PCP 3 (\i(x),j3i =l,a 


= 0.025) 


0.667 1.000 


0.961 


1.000 


1.000 


PCP 3 {X l {x),p 1 = l,a 


= 0.05) 


0.208 1.000 


0.625 


0.947 


0.927 


PCP 3 {X l {x),p 1 = l,a 


= 0.1) 


0.091 0.968 


0.338 


0.570 


0.560 


PCP 3 (X 1 (x),p 1 =2,a 


= 0.025) 


0.187 0.789 


0.458 


0.495 


0.411 


PCP 3 (X 1 (x),p 1 =2,a 


= 0.05) 


0.077 0.216 


0.129 


0.129 


0.112 


PCP 3 (X 1 (x), f 3 1 =2,a 


= 0.1) 


0.042 0.089 


0.079 


0.073 


0.069 


PCP 3 (X 2 (x),[3 2 = 3,a 


= 0.025) 


0.226 1.000 


1.000 


1.000 


1.000 


PCP 3 (X 2 (x),f3 2 = 3,a 


= 0.05) 


0.145 0.976 


1.000 


1.000 


1.000 


PCP 3 (X 2 (x),f3 2 = 3,a 


= 0.1) 


0.089 0.776 


1.000 


1.000 


1.000 


PCP 3 {\ 2 (x),f3 2 =5,a 


= 0.025) 


0.826 1.000 


1.000 


1.000 


1.000 


PCP 3 {\ 2 {x),p 2 =5,a 


= 0.05) 


0.476 1.000 


1.000 


1.000 


1.000 


PCP 3 {X 2 (x),r3 2 =5,a 


= 0.1) 


0.385 1.000 


1.000 


1.000 


1.000 



Table 2: Power of the integral deviation measures D according to the isotropic (I), border 
(B), modified border (MB), translation (T) edge correction factor and without edge 
correction (N). 



4 Influence of the intensity estimate 

The estimators of the second-order characteristics defined in Section 12.51 are approximatively 
unbiased. However, they depend on the first-order intensity function A(x). As X(x) is not 
known in practice, it must be replaced by an estimate, A(x), in Equatio ns (Bl) and Q . 
There are s everal wa ys in estimating the i ntens ity function; see for example ICressid (J1993I ). 
Illian et all ( 20081 ) or IStoyan and Stoyanl ( 19941 ) for a review. The two mains are using a 



parametric model, often used when the point pattern suggest a theoretical form of A(x), or 
using kernel methods. In this section, we illustrate by simulations the influence of X(x) on 
the estimation of the pair correlation function. 
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4.1 Estimation 

In the following, we assume the spatio-temporal separability of X(x). Thus, whatever the 
method used, the estimation of the intensity is treated completely separately for space and 
time. 

4.1.1 Parametric estimation 

We consider an exponential model for the tempor al intensity func t ion and a log linear functio n 
of spatial coordinates for the spatial intensity (JLiu et all . 120071 : IStoyan and Stoyanl . Il994l ): 
A(s; (3) = exp(/? T Z(s)), where Z(s) is a vector of polynomials of coordinates s x ,s y of point 



location s. Polynomials are of order one and two: 

Model 1: A(s; j3) = exp(/3 + /3is x + /3 2 s y ). 

Model 2: A(s;/3) = exp(/3 + fcs x + <3 2 s y + fas 2 x + PiS x s y + /3 5 aJ). 

Parameters are estimated using likelihood methods. The log-likelihood for (3 is 

n ,. 

logL(/3,s) = ^log(A(s 4 ;/3))- / A(u;/3)du. 



We use the Berman- Turner algorithm ( Berman and Turner! . 1 19921 ) for finding the maximum 
likelihood estimate of /3. 



4.1.2 Kernel estimation 



The temporal intensity function is estimate d using a gaussia n kernel estimator, with band- 
width 0.9n -1 ' 5 min{<r, (Q3 — Qi)/1.34}, see ISilvermanl (119861 ). To estimate the spatial inten- 



sity function we use a quartic kernel estimator. The estimator takes the form 

kh(s ~ Si) 



k*) = E 



csM s i) 



where kh{s) = h 2 k(s/h) and cs t h(si) = §s ^h(s-Si) ds is the edge correction factor defined in 



error 



Diggle (119851) ensuring that f„ M s) ds = n. The bandwidth, h, minimises the mean squared 



Berman and Diggld . Il989l ). 



4.2 Simulations 

Simulations are based on the 1000 realisations of the inhomogeneous Poisson processes 
IPP(\(x),f3) and of the Poisson cluster processes PCP^{X{x),f3, a) described in Section [3721 
According to the results of Section 3, we use the translation method to correct edge effects 
and consider the pair correlation function. For each simulation, we compute the pair correla- 
tion function according to the different estimates of the intensity function. In the following, 
the estimator is indexed by A when using the true intensity, by A p when using a parametric 
estimate and by A& when using a kernel estimate of the intensity function. 
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4.3 Performance of g relative to A 

We compute the pointwise bias and mean squared error of the pair correlation function esti- 
mator according to the different estimates of the intensity function. For the inhomogeneous 
Poisson process IPP(\(x),/3), g~\ and g^ have similar MSE and negative bias, both increas- 
ing with u and v. Bias and MSE of gT are greater than for the two others, maybe due to 
over-smoothing, and much more greater for small values of u and v and as the inhomogeneity 
is stronger. The bias is positive for small u and v and becomes negative for large u and 
v. This is illustrated in the first two rows of Figure [U The last two rows of Figure 0] show 
results for the Poisson cluster processes. Model 2 implies very high bias at different ranges of 
spatial and temporal distances, thus leading to very high values of the corresponding MSE. 
For PCPs(\i(x),^i,a) the bias and MSE of g\ and g^ (Model 1) have similar orders of 
magnitude, except for small u, where they are greater when using parametric estimates of 
the intensity. This is not the case for PCP^{\2{x),f32^)- The estimator <fT is biased for 
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Figure 4: Bias of the pair correlation function estimated by using A (first column), Afc (second 
column), A pi (third column) and A P2 (fourth column). Spatial distances u are in 
abscissa, while temporal distances v are in grey with the same values than u (the 
darker, the greater). 

the Poisson cluster processes. This is mainly due to the fact that we estimated non paramet- 
rically both the intensity function and the pair correlation function from the same observed 
point pattern. Without any other assumption or information about the nature of the under- 
lying point process, we cannot in general make an empirical distinction between first-order 
and second-order effects and thus cannot distinguish the contributions due to spatial inho- 
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mogeneity or spatial interacting phenomena. Here, the bandwidths h selected for spatial 
kernel smoothing (see Section I4.1.2P are closed to the spatial dispersion a. This leads to 
over-smoothed intensity functions and thus to severely biased values of the pair correlation 
function. 



4.4 Clustering detection relative to A 

We analyse the detection of spatio-temporal clustering as in Section 13.51 The probability of 
detecting clustering is computed for g~\, g^ and g^ from realisations of the Poisson cluster 
processes PC P%(\(x) , j3 , a) . One way to get around the difficulty of estimating both the 
first and second-order characteristics non-parametrically consists in ass uming that firs t -order 
effects oper a te at larger spatial scale than the second-order effects (jBaddeley et all . l2000l : 
Diggle et all . 120071 ). Thus, g^ is also evaluated for fixed spatial bandwidths, h = 0.2 and 
h = 0.4. For the processes PCP 3 (Ai(x),/3i = I, a) and PCP 3 (X 2 (x), fa = 3, a), Figure E 
shows the probabilities of detecting clustering when using the true intensity (left), the kernel 
estimation with bandwidth h = 0.4 (middle) and the parametric model 1 (right). It appears 
that clustering is well detected for g\ and g^ . We get a lower detection rate for g^ . 
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Figure 5: Probabilities of detecting clustering using the true intensity function (first column), 
the non-parametric estimate with bandwidth h = 0.4 (second column) and the 
parametric estimate from Model 1 (third column) for the Poisson cluster processes 
PCP 3 (\(x),l3,a) with a = 0.025. 

Finally, we calculate for each estimator the overall statistic D and the related probabilities 
of detecting clustering. Because D accumulates deviations of g from the theoretical g and 
because the bias of g^ and g^ may have high values, to avoid misleading results for these two 

estimators we rather consider the analog of D defined in Equation ([S]) . Table shows these 
probabilities for various scales of spatial clustering a and inhomogeneity. In the simulated 
Poisson cluster processes, the first-order effects operate at larger spatial scale than do the 
second-order effects, so we can see that considering this assumption clearly improves the 
results and that spatio-temporal clustering can now be detected. The estimators g^ and g^ 
are sensitive to a and (3; their power decrease as a and /3 increase, except when the bandwidth 
h is used. Indeed, it increases with a, thus leading to less biased estimates. Conversely, the 
performances for a fixed bandwidth decrease with a because it becomes then more difficult 
to distinguish first- and second-order effects. 
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h 
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h = 0A 


Model 1 


PCP 3 (X x (x),^ 
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1.000 


0.000 


0.999 


1.000 


0.999 


PCPMx),^ 
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1.000 
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PCP 3 (X 2 (x),f3 2 


= 3,C7 = 0.1) 


1.000 
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0.596 


0.994 


0.288 


PCP 3 (X 2 (x),(3 2 


= 5,(7 = 0.025) 


1.000 


0.003 


0.481 


0.801 


0.972 


PCP 3 (X 2 (x),p 2 


= 5,(7 = 0.05) 


1.000 


0.126 


0.191 


0.334 


0.604 


PCP 3 (X 2 (x),r3 2 


= 5,cr = 0.1) 


1.000 


0.268 


0.096 


0.107 


0.188 



Table 3: Power of the integral deviation measures D according to the method used to esti- 
mate the intensity function: kernel estimation with auto-selected bandwidth h or 
fixed bandwidth h; parametric estimation (Model 1). 

5 Discussion 



The spatio-temporal inhomogeneous ET-function and pair correlation function describe second- 
order characteristics of point processes. They can be used as an exploratory tools, e . g. for 
testing spatio-temporal clus tering or spatio-temporal interaction ( Diggle and Gabriell . I201CJ ; 
Moller and Ghorbanil . l2012l ). Their non-parametric estimation depends on an edge correc- 
tion factor and assumes that the intensity function is known. We extended classical spatial 
edge correction factors to the spatio-temporal setting and explored the influence of edge cor- 
rection methods and intensity estimation on the performance of these estimators. Results 
on edge effects indicate that the border method performs well when the point pattern may 
be assumed homogeneous, otherwise, as soon as the point pattern tends to be inhomoge- 
neous/clustered/anisotropic the translation method must be preferred. Some preliminary 
tests for isotropy and Poisson assumption shoul d be done before estimating second-order 
characteris t ics; se e for instance iGuan et all ( 200a ) for testing isotropy in t he spat i al cas e or 
Illian et all ( 20081 ) for some discussion on the detection of anisotropy, and iDiggld ( 20031 ) for 
testing the Poisson assumption. In the case of anisotropy, one then can apply adapted ver- 
sion of the iv-function and pair correlation (jMoller and Toftagerl . |201_2]; IStoyan and Stoyanl . 
19941 ) . Note that estimators without edge correction are severely biased for values of the 
spatial and temporal distances larger than the ones used in this study. This is particularly 
so for the estimator of the STIK-function. Results related to the intensity estimation show 
that the performance of the pair correlation function can be severely altered by the inten- 
sity estimate. This can be explained by over-parametrisation or over-fitting in the case of a 
parametric estimation of the intensity function or by the incapacity of distinguish first- and 
second-order effects from a single realisation of the point process in the case of a kernel-based 
estimation. Non-parametric estimation of first- and second-order characteristics requires the 
assumption that first- and second-order effects operate at different scales. So, an important 
question is what to do if first- and second-order effects operate at the same scale? It would 
be interesting to understand where does the confusion between the effects come from. 
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Appendix 1 

The proof of the unbiasedness property of the STIK-fu nction estimator given in Equation ([3]) 
is based on the following Proposition, see for example iMoller and Waagepetersenl (120031 ). 



Proposition 1 Suppose that a point process X defined on W C M. d has a second-order 
intensity X 2 - For a function f : M d x M rf — > [0, 00), 



E 



+ 



x,y£X w 



f(x,y)X 2 (x,y)dxdy. 



(9) 



For VT = 5xTgIR 2 x1R and x = (s x ,t x ), y = (s y , t y ), we have from Equation © 



* 1 1 

k(u, v) = y^ — — - . 



x,y<=X w 



w(x,y) X(x)X(y) 



~Sy\\<U ; \t X ~ty \<V\- 



Thus, for h = x—y = (h s , h t ) and f(x, y) = x f ^ N - j—^ ( v ^ we have from Equation ([9]) 



E 



K(u,v) 



1 A 2 (z,2/) 
w(x,y) X(x)X(y) 



w(x,y)X(x)X(y) 

-s y \\<u ; \t x -t y \<v} 1 w(x)lw(y)dxdy 



1 



w(x,x + h) 



lw(x)lw( x + h) dx g(h)ln 



h s \\<u ; |/it|<u} 



dh 



(*) 



For the isotropic edge correction method, denoting wt(x, y) = 1 if [t x — |t x — t y \; t x + |t x — t y \]n 
T 7^ and 1/2 otherwise, and |<%(s,r)| the circumference of the disc of center s and radius 
r, 



(*) 



1 



|56(Sa:, /l s )| 1sxt(Sx, t x )l S xT(s x + h s ,t x + h t ) 



\S xT\ J J \db(s x ,h s )nS\ 
f \db(s x ,h s 



w t (x,x + h) 



ds x dt x 



S xT\ J \db{s x ,h s )r\S\ 



ls(s x )ls(s x + h s )ds x f 



l T (t x )l T (t x + h t ) 



w t (x,x + h) 



dU 



1. 



For the modified border method, 
M = K Vl^l / 1 {Seu*T ev }{x)lw{x + /»)l{||h,||<« ; |h,|< t ,}dx 



~ T I / 1 {S eu y.T ev }{x)^W eh {x)l{\\ ha \\< u ■ \ ht \<v}dx 



1 



Pea x Tq v 



L {(5e«nSeOx( T e..nTefc t )}( a; ) 1 {|IM<« ; |fc*l<«> 



d.r 



xT I / ^Se^xTe,}^)^ = 1. 

u x J 0t>| J 
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For the translation method, 
(*) = 



1 



sns hs \ x\Tnr ht 
1 



\sns hs \ x\TnT ht \ 



1sxt(sx, t x )lsxT(s x + h s ,t x + h t ) ds x dt x 
1 sns^ h3 {sx)^-Tf\T^ ht (t x ) ds x dt x = 1- 



Thus, 



E 



K(u,v) 



9{h)l{\\h s \\<u ; \h t \<v}dh 

g(h s ,h t ) dh s dh t = K(u, v) 



v Jb(0,u) 



Under the assumption of isotropy, we obtain Equation ([2]) by using cylindrical coordinates of 
9(h). 
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Appendix 2 

Figure [6] shows the mean squared error of the pair correlation function estimated when using 
the isotropic (I), border (B), modified-border (MB), translation (T), none (N) edge correction 
methods. Spatial distances u are in abscissa, while temporal distances v are in grey: the 
darker, the greater, with the same range of values than u. 
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Figure 6: Mean squared error of the pair correlation function from IPP(\i(x), /3± = 1) (first 
row), PCPi{<7 = 0.05) (second row), PCP 3 (Ai(x),/3i = l,a = 0.05) (third row), 
PCP 3 (\ 2 (x),fi2 = 3,cr = 0.05) (fourth row), PCP 2 (a = 0.2, C 2 = 0.0625) (fifth 
row) and PCP 2 {(J = 0.2, ( 2 = 0.5) (sixth row). 
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Point process 



HPP(X) 
IPP(Xi(x),$i = l) 
IPP(\ 1 (x),p 1 = 2) 
JPP(A 2 (x),/3 2 = 3) 
/PP(A 2 (x),ft, = 5) 
PCPi(o = 0.025) 
PCPl(cr = 0.05) 
PCPi(cr = 0.1) 
PCPi(cr = 0.15) 
PCPi(cr = 0.2) 
PCP 3 (Xi(x 
PCP 3 (Xi(x 
PCP 3 (Xi(x 
PCP 3 (Ai(x 
PCP 3 (Ai' 
PCP 3 (Ai 
PCP 3 (Ai 
PCPsCAi 
PCP 3 (Ai 
PCP 3 (Ai 
PCP 3 (A 2 
PCP 3 (A 2 
PCP 3 (A 2 
PCP 3 (A 2 
PCP 3 (A 2 
PCP 3 (A 2 
PCP 3 (A 2 , 
PCP 3 (A 2 (x 
PCP 3 (A 2 (z 
PCP 3 (A 2 (x 



pi = l,o- 
Pl = l,o-- 

Pl = l,<T- 

/3l = l,<j: 
Pi = 1,<T: 
Pl = 2,<T ■ 

$i = 2,<t-. 

/3l=2,0-: 

Pl = 2,o--. 
$1 = 2,0--. 
/3 2 = 3,o-- 
$2 = 3, a ■- 
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$2 = 3,0--. 
$2 = 3,<J- 
$2 = 5,o--. 
$2 = 5,o--. 
$2 = 5,o--. 

$2 = 5,0-: 
$2 = 5,0-: 



0.025) 

0.05) 

0.1) 

0.15) 

0.2) 

0.025) 

0.05) 

0.1) 

0.15) 

0.2) 

0.025) 

0.05) 

0.1) 

0.15) 

0.2) 

0.025) 

0.05) 

0.1) 

0.15) 

0.2) 



1 



Edge 
B 



correction 
MB 



factor 



N 



5.56 
5.30 
2.95 
5.21 
3.74 
0.41 
0.20 
0.08 
0.18 
0.30 
2.02 
1.88 
1.37 
1.69 
1.55 
3.36 
1.08 
0.44 
1.09 
0.67 
0.77 
0.52 
0.35 
0.25 
0.24 
0.64 
0.40 
0.17 
0.14 
0.10 



100.00 

100.00 

100.00 

100.00 

100.00 

1.15 

0.67 

1.13 

1.63 

4.85 

20.46 

21.61 

8.60 

19.57 

26.00 

71.25 

60.08 

65.55 

64.23 

80.70 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 



1.42 

2.13 

4.35 

11.37 

9.23 

0.16 

0.15 

0.26 

0.34 

0.99 

11.94 

15.47 

6.28 

22.76 

18.06 

100.00 

100.00 

100.00 

100.00 

100.00 

1.29 

0.75 

0.65 

0.42 

0.44 

1.75 

1.96 

0.79 

0.69 

0.93 



13.10 

7.20 

6.83 

54.57 

14.84 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

61.53 

19.78 

7.36 

31.25 

10.44 

18.44 

14.38 

18.05 

12.85 

11.02 

5.73 

5.79 

4.13 

3.64 

2.38 



22.09 

12.73 

12.86 

64.74 

27.22 

13.92 

9.26 

9.62 

25.01 

57.11 

37.28 

30.14 

21.11 

26.05 

26.71 

20.90 

7.13 

2.54 

8.06 

3.18 

8.93 

6.17 

6.78 

4.84 

4.13 

3.09 

2.70 

1.76 

1.52 

0.99 



Table 4: Relative efficiency of the integral deviation measures D based on the STIK-function 
according to the isotropic (I), border (B), modified border (MB), translation (T) 
edge correction factor and without edge correction (N). 
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Edge correction factor 




Point process 


I 


B 


MB 


T 


N 


PPP(A) 


0.79 


100.00 


10.05 


15.60 


23.21 


IPP(\l(x), fa = 1) 


4.22 


100.00 


24.88 


18.11 


31.56 


IPP(\ 1 (x),/3 1 = 2) 


5.68 


47.31 


60.84 


100.00 


48.77 


7PP(A 2 (x),/3 2 = 3) 


0.01 


26.07 


23.90 


100.00 


59.21 


IPP(X 2 (x), /9 2 = 5) 


0.03 


44.02 


41.72 


100.00 


86.46 


PCPi(o = 0.025) 


3.80 


9.09 


17.97 


100.00 


28.73 


PCPi(o- = 0.05) 


1.75 


7.71 


12.54 


100.00 


27.70 


PCPl(cr = 0.1) 


0.32 


4.66 


6.49 


100.00 


22.43 


PCPilo = 0.15) 


0.36 


4.96 


3.22 


100.00 


33.64 


PCPi(<t = 0.2) 


0.39 


6.05 


5.96 


100.00 


42.57 


PCP 2 (a = 0.025, C 2 = 0.0625) 


17.26 


28.64 


44.93 


100.00 


55.38 


PCP 2 \o = 0.05, C 2 = 0.0625) 


10.76 


24.72 


43.67 


100.00 


51.26 


PCP 2 (o = 0.1, C 2 = 0.0625) 


5.42 


18.40 


24.01 


100.00 


42.46 


PCP 2 (o- = 0.15, C 2 = 0.0625) 


2.59 


12.72 


19.71 


100.00 


36.33 


PCP 2 (o = 0.2, C 2 = 0.0625) 


1.56 


10.67 


20.48 


100.00 


32.35 


PCP 2 (ct = 0.025, C 2 = 0.25) 


23.34 


32.08 


50.91 


100.00 


58.44 


PCP 2 (cr = 0.05,C 2 = 0.25) 


14.71 


27.28 


44.38 


100.00 


53.93 


PCP 2 (o = 0.1,C 2 =0.25) 


6.21 


20.36 


27.30 


100.00 


45.17 


PCP 2 (o- = 0.15,C 2 =0.25) 


4.24 


14.07 


18.11 


100.00 


39.83 


PCP 2 (o- = 0.2,C 2 =0.25) 


2.15 


12.97 


19.72 


100.00 


35.26 


PCP 2 (a = 0.025, C 2 = 0.5) 


26.39 


34.33 


52.15 


100.00 


60.31 


PCP 2 (cr = 0.05,C 2 =0.5) 


18.39 


31.36 


50.26 


100.00 


57.38 


PCP 2 (o = 0.1, C 2 =0.5) 


13.32 


26.29 


35.91 


100.00 


51.86 


PCP 2 [o = 0.15,C 2 =0.5) 


5.50 


18.39 


27.61 


100.00 


46.79 


PCP 2 (o- = 0.2, C 2 =0.5) 


4.44 


16.93 


19.81 


100.00 


42.37 


PCP 2 (ct = 0.025, C 2 = 0.75) 


26.53 


36.53 


54.16 


100.00 


61.38 


PCP 2 (cr = 0.05,C 2 = 0.75) 


23.29 


33.84 


54.01 


100.00 


59.30 


PCP 2 (o = 0.1,C 2 =0.75) 


17.11 


29.93 


41.53 


100.00 


54.43 


PCP 2 {o = 0.15,C 2 =0.75) 


6.87 


21.24 


36.98 


100.00 


50.36 


PCP 2 {o = 0.2,C 2 =0.75) 


9.50 


20.00 


24.96 


100.00 


47.72 


PCP 2 (o- = 0.025, C 2 = 1) 


32.54 


36.99 


53.98 


100.00 


61.74 


PCP 2 (cr = 0.05, C 2 = 1) 


24.98 


34.67 


53.14 


100.00 


59.70 


PCP 2 [o- = 0.1, c 2 = 1) 


3.52 


31.16 


44.69 


100.00 


56.02 


PCP 2 (cr = 0.15, C 2 = 1) 


8.60 


23.10 


42.20 


100.00 


52.47 


PCP 2 (cr = 0.2,C 2 = 1) 


10.22 


22.72 


30.69 


100.00 


49.94 


PCP 3 (AiO),/3i = l,o- = 0.025) 


6.70 


7.29 


74.58 


100.00 


46.59 


PCP 3 (Ai(x),^i = l,o- = 0.05) 


3.36 


6.03 


47.78 


100.00 


33.69 


?CP 3 (Ai(4ft = 1,(7 = 0.1) 


2.71 


5.79 


15.02 


100.00 


29.54 


PCP 3 (Ai(x),/3i = l,cr = 0.15) 


0.36 


6.80 


19.07 


100.00 


30.72 


PCP 3 (Ai(x),/3i = l,cr = 0.2) 


1.99 


7.53 


23.50 


100.00 


26.82 


PCP 3 (Ai(x),/3i =2,(7 = 0.025) 


12.90 


16.90 


59.70 


100.00 


45.67 


PCP 3 (Ai(i),/3i = 2,o- = 0.05) 


12.63 


41.33 


76.29 


100.00 


47.71 


PCP 3 (Ai(z),/3i = 2,cr = 0.1) 


11.78 


97.97 


69.36 


100.00 


43.24 


PCP 3 (Ai(z),/3i =2,cr = 0.15) 


5.94 


47.25 


75.73 


100.00 


38.38 


PCP 3 (Ai(x),^i = 2,(7 = 0.2) 


4.76 


100.00 


31.79 


62.56 


21.84 


PCP 3 (A 2 (z),/3 2 = 3,cr = 0.025) 


0.24 


100.00 


4.82 


9.05 


5.62 


-PCP 3 (A 2 (z),,3 2 = 3,cr = 0.05) 


0.28 


100.00 


2.39 


6.36 


3.67 


PCP 3 (A 2 (x),/3 2 = 3,(7 = 0.1) 


0.06 


100.00 


1.67 


6.87 


3.43 


-PCP 3 (A 2 (x),/3 2 = 3,(7 = 0.15) 


0.05 


100.00 


0.90 


3.88 


1.93 


PCP 3 (A 2 (x),& = 3,(7 = 0.2) 


0.06 


100.00 


0.74 


3.13 


1.56 


PCP 3 (X 2 (x),l3 2 = 5,cr = 0.025) 


0.29 


100.00 


2.56 


4.48 


3.24 


PCP 3 (A 2 (x),/3 2 = 5,(7 = 0.05) 


0.23 


100.00 


2.13 


3.85 


2.54 


PCP 3 (A 2 (x),& = 5,(7 = 0.1) 


0.10 


100.00 


1.19 


3.24 


1.94 


PCP 3 (A 2 (z),/3 2 = 5,cr = 0.15) 


0.04 


100.00 


0.84 


2.50 


1.46 


PCP 3 (A 2 (x),ft = 5,(7 = 0.2) 


0.10 


100.00 


0.87 


2.18 


1.27 



Table 5: Relative efficiency of the integral deviation measures D based on the pair corre- 
lation function according to the isotropic (I), border (B), modified border (MB), 
translation (T) edge correction factor and without edge correction (N). 
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